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SUMMARY 


An investigation has been made of the use of on-off reaction jets 
for precision attitude control of a satellite. Since a symmetrical 
vehicle is assumed, only single-axis control needs to be considered. 

The responses to initial disturbances and also limit-cycle character- 
istics for several systems have been evaluated. Calculated results 
indicate that realistic values of settling time and fuel consumption 
for the example considered can be obtained. 

The performance of a given system depends on the characteristics 
of the error detector used. In cases where the detector output was 
saturated for a relatively low error input, the settling time deterio- 
rated when a lead network was used to provide damping. This deterio- 
ration could be eliminated if a separate rate signal to produce vehicle 
rate limiting were available. As an alternate approach, two systems 
were investigated which used a timed sequence of torques and could 
operate with a detector output of very small linear range. Although 
the performance of these systems was poorer than that of the lead net- 
work system without detector saturation, the performance was better 
than that of the lead network system with low values of detector satura- 
tion. The effects on limit-cycle characteristics of hysteresis, lead 
network constants, dead zone, and thrust time delays were also investi- 
gated. 


INTRODUCTION 


As more sophisticated satellite payloads become feasible, more pre- 
cise attitude control systems generally become necessary. An example of 
this trend is the orbiting astronomical observatory (ref. l). The atti- 
tude control requirements for this project include space pointing accu- 
racies as high as 0.1 second of arc as well as relatively high angular 
rate capability in order to slue from one target star to another. In 
addition, the target may be occulted by the earth for part of the orbit, 
and the vehicle must be reoriented when the target reappears. 
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General aspects of the attitude control, problem were discussed in 
references 2 and 3* The use of inertia wheels for precise attitude con- 
trol is discussed in reference 4. An alternate approach using on-off 
reaction jets, which is attractive because of its simplicity of oper- 
ation, will be discussed in this report. References 5 and 6 discussed 
principally the limit-cycle characteristics of on-off jet control sys- 
tems for cases with pointing accuracy requirements less severe than 
those considered here. Investigations have also been made of many of 
the characteristics of on-off control systems which use a linear combi- 
nation of state variables to form a switching control function (ref. r j) 
or which use a nonlinear combination of state variables to produce a 
system with minimum response time to a giver input or initial condition 
(ref. 8). 

For the orbiting astronomical observatory, both a coarse and fine 
control mode are used because of the large dynamic range required. 

This report will be concerned with the characteristics of a fine con- 
^'^’ < - ) l mode; that is, it will be assumed that a coarse control mode has 
alined the vehicle to within a set of prescribed error and error rate 
values. The use of microthrust jets to correct these initial errors, 
and the resulting low-frequency limit-cycle motion about a prescribed 
dead zone will be investigated. Since the vehicle is inertially 
balanced as closely as possible, and since tie reaction jets do not 
introduce cross-coupling torques, coupling effects between control axes 
w 411 be negligibly small. Hence the performance can be analyzed on a 
single degree of freedom basis. The influence of external torques due 
to solar, gravity, magnetic, and atmospheric effects must also be con- 
sidered. These disturbances vary at rates associated with the orbital 
period. Since system settling times of the srder of 3 minutes are de- 
sired, the external torques may be treated a ; constants. 

An additional important effect which wi.l be investigated is that 
of detector saturation. This effect may be encountered for any system 
for which a large dynamic range is required. While an error signal 
varying with amplitude over the entire control range is certainly desirable 
even though multiple detector ranges would be required, certain geometric 
and sensitivity limitations may make a continuously varying signal 
difficult to obtain. Two systems which use a timed torque sequence and 
can operate in the presence of detector saturation will be investigated. 
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SYMBOLS 


magnitude of angular acceleration applied by reaction jets or 
average value if positive and negative values are unequal 

magnitude of negative acceleration applied by reaction jets 

magnitude of positive acceleration applied by reaction jets 

acceleration produced by external torque 

force applied by reaction jets 

relay total hysteresis 

distance from reaction jets to center of gravity 
moment of inertia of satellite vehicle 
specific impulse 

gain of differentiating network used with linear switch-curve 
system 

transform variable 
time^ sec 

time required to complete one cycle for timer or rate systems 

delay time approximating thrust decay 

delay time approximating thrust rise 

total settling time for a given initial condition 

average total settling time for a range of initial conditions 

total settling time for an initial angular rate 

time constant for differentiating network 

fuel weight consumed during a given transient motion 

average fuel consumption rate during limit-cycle motion 

average normalized fuel consumption rate during limit-cycle 
motion, per sec 
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sgn( ) 


normalized, fuel weight required for one cycle of tinier or rate 
systems 

W I l 

normalized fuel weight for a giver transient motion, Ip- 

average normalized fuel weight for a range of initial con- 
ditions 

angular displacement 

do 

angular velocity, — 

CL U 

maximum angular deviation during coast portion of limit cycle 
with no external torque 

desired control accuracy of attitude angle 

control function dead zone of linear switch-curve system 

desired control accuracy of angular rate 

saturation value of error detector 

angular- rate limit 

measured value of 9 at 6 = 0 for rate system 

increment in 0 produced "by jet Ihrust during limit-cycle 
motion 

torque-on time for motions near ti e limit-cycle condition 
torque-on time for limit cycle 

k )i 
( ) 


Subscript 


o 


initial value 



ANALYSIS 


The characteristics of several systems which use low- thrust on-off 
jets for the attitude control of a satellite, such as the orbiting astro 
norriical observatory, will be developed in this section. As stated in 
the introduction, only a single degree of freedom will be considered. 

The system will be assumed to be alined initially within a given value 
of error and error rate. Equations will be developed for determining 
settling times and fuel consumption for systems with different forms of 
switching logic to control the jet action. The effects of detector 
saturation on performance will also be included. In addition, equations 
for limit-cycle characteristics and the corresponding fuel consumption 
will also be determined. Since transient performance is of primary 
interest, it will be convenient to illustrate vehicle response through 
use of the phase plane (e.g., ref. 9)* 

For an application, such as the orbiting astronomical observatory, 
the initial errors are large in comparison with the ultimate tracking 
accuracy desired even if separate coarse and fine control modes are 
used. For instance, if an initial error of 2 minutes of arc for a fine 
control mode is assumed, with a desired tracking accuracy of 0.1 second 
of arc, the initial error would be over 1000 times the desired value. 
Because of the large dynamic range being considered, the analysis can be 
separated into two parts: (l) the settling time required for the sys- 

tem to recover from an initial error to within a desired value, and (2) 
the transition into a limit-cycle motion. For the astronomical observa- 
tory program, every effort is made to reduce the magnitude of the ex- 
ternal torques because of the high pointing accuracy required. Initial 
estimates were that by careful balance of inertial, magnetic, solar re- 
flecting, and aerodynamic properties, average torques could be reduced 
to the order of 100 dyne-cm for a vehicle with moments of inertia of 
about 10 lo gm cm 2 (ref. l). Although this value of torque is small in 
comparison with later estimates due to subsequent changes in the vehicle 
it will be adequate for use in the example calculations which will be 
shown in this report. Since the applied torque to be considered will be 
about 100 times larger than the external torque, the latter will affect 
only the limit-cycle motion. Similarly, other properties such as dead 
zones and time delays significantly influence only the limit-cycle 
characteristics. However, it should be noted that calculations of time 
to settle to within a given error could be affected If this error were 
sufficiently close to errors encountered during limit-cycle motion. 

Since external spring and damping torques acting on the controlled 
vehicle can be neglected, the resulting vehicle motion which results 
from the control system applied torques can be given approximately in 
analytical form. An analog computer was used mainly as a check on the 
calculated results from the equations which will be developed. 



Recovery From Initial Errors 


Common characteristics of the systems to be considered will be 
discussed first. A block diagram is shown in figure 1. The angular 
error is measured by an error detector with assumed perfect dynamic 
characteristics. The switching function, for which several types will 
be compared, determines the direction of the torque to be applied to 
the vehicle. The function actuates the appropriate relay which controls 
an on-off type valve that applies a torque in the proper direction. The 
resulting equation of motion is linear during the interval the applied 
torque is constant. When the applied torqoe switches sign, continuity 
in 0 and 0 is maintained so that initial conditions for the suceed- 
Ing motion are the same as the end conditions for the previous motion. 
Since the vehicle damping and spring forces can be neglected, and since 
the applied torque is much greater than the assumed disturbance torques, 
the applied torque FZ, produces a constant acceleration 


•• FZ 

0 = ±a = =j- 

(1) 

The resulting vehicle motions become 


6 = 9 0 — at 

(2) 

0 = 0 O + 0 o t ± (i/;?)at 2 

(3) 


The zero subscripts refer to conditions existing when the torque is 
initially applied. 

If no torque is applied, the vehicle coasts at a constant angular 

rate 


0 = 0o 

(*0 

0 = 0o + 0 o t 

(5) 


The zero subscripts for this case refer to initial conditions at the be- 
ginning of the coast period . 

The switching times are determined by functions of the state of the 
system which will be described by the variables 0 and 6 . A convenient 
way of representing these motions is through use of the phase plane 
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(ref. 9). 'When time is eliminated between equations (2) and (3)> the 
motions in the phase plane become (l) for the constant torque case 


e = e 0 


iol . ©L 

+2a +2a 


( 6 ) 


and (2) for the coast case 

0 = 0 O 


Complete motion histories can be obtained by piecing together these 
trajectory segments. 

To provide a quantitative comparison of different systems, two 
parameters are of interest: (l) settling time, that is, the time re- 
quired for the system to recover and remain within a certain desired 
amplitude and (2) the fuel consumption during this motion. The total 
settling time required, t , for a given motion is determined from a 
series of trajectory segments as 


J 

-E 


A0 

IV 

A0 

a 

+ E 

6 


k=i 


(7) 


where J is the number of torque-on trajectory segments, K is the number 
of torque-off trajectory segments, and A( ) denotes the change in a 
variable over a trajectory segment. The weight of the fuel consumed 
during the motion can be obtained (assuming a constant specific impulse 
for all operating conditions) by the sum of the changes in angular rate 
of the torque- on trajectory segments. 


where 


IW t 

W = - 

til 

sp 


w, = 


£ m 


J = 1 


( 8 ) 


Settling time and fuel consumption for several systems will be 
compared. Since a nonlinear system may be tuned to a particular set of 
initial conditions, a range of sets of initial conditions, each suf- 
ficient to define the initial state of the system, must be investigated. 
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For this report initial values of angle and angular rate errors will he 
used. 1 For each system, an average value of settling time and fuel 
consumption will be determined, based on the :’ange of initial conditions 
to be considered. The selected conditions should be as realistic as 
possible to provide a valid comparison between the systems. A simpli- 
fied set of initial conditions, which are distributed along a line and 
normalized to a desired angular tracking accuracy 0^, that was used 
for this report is shown in sketch (a). 
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Sketch (a) 

This distribution is representative of that which could be introduced by 
a coarse-control on-off system for which the switch curve is linear. In 
addition, the probability of occurrence of al] combinations of Q 0 / e & 
and ^o/^d was assumed to be equal. After a direction of integration, 
dx, is defined along the line of initial conditions, and the probability 
of occurrence of the initial conditions is set equal to unity so that 
the length of the line in the direction x is unity, the previous 
equation reduces to 



t 


t 



(9) 


1 If a joint probability distribution of initial values of angle and 
angular rate is known, p(0oj>0o)> then an average value of settling time 
can be determined by the following integratior. over a surface, S, in the 
Qo> 0o plane, 

t t = / s /p( 6o> e O )t(0 o ,S o ) d0 c 

The average settling times can then be compared for several systems. 


< -3- cr\ oo 
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The integration was performed numerically along the line describing the 
initial conditions shown in sketch (a). The average values of settling 
time, t , were obtained by evaluating equation ( 9 ) hy use of Simpson’s 
rule from calculated values of t^ at equally spaced increments of 
0o,/0d and Oo/^d * The average values of normalized fuel weight* 

Wf/Sd* were calculated by a similar numerical integration. 

The equations used to determine system settling time and fuel con- 
sumption during recovery from an initial error will be developed next. 

Linear switch curve system. - For this system* the switching logic 
used is of the following form: 








+ 


K 0 S 





Sketch (b) 


The rate signal is assumed to be formed by a perfect lead network. The 
effect of the time constant of the network was found to be important 
only during limit-cycle motion and, hence* was neglected. The resulting 
motion is then determined by the equation 


9 = -a sgn(0 + K“0) 

w 

Thus* the torque changes sign whenever the quantity 0 + K ^0 changes 
sign. For convenience* a unit gain is associated with 6 . 

To obtain the total settling time* the motions shown in sketch (c) 
are separated into three parts: 

(1) The time required to first reach the switch curve from a given 
set of initial conditions* 

(2) The time required during which the trajectory overshoots the 
switch curve when it is intersected* and 

( 3 ) The time during which the motion subsides as commanded by the 
switch curve until the desired attitude angle is reached. 
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Sketch (c) 

Equations for total settling time and fuel consumption may be obtained 
from equations for the trajectory segments, the switch- curve equation, 
and relations describing the time required and fuel consumption for 
each trajectory segment. A derivation of the equations is given in the 
appendix. The equations to be given in this report for settling times 
and fuel consumption will be valid for a range of initial conditions at 
least within the first quadrant of the phase plane. This is sufficient 
to evaluate the initial conditions shown in sketch (a). While equations 
could be written for more general initial conditions, they tend to be- 
come rather cumbersome. 

The time, ti, required to reach the switch curve from initial con- 
ditions G 0 , 6 o 



The total time required to settle to a given amplitude 0 can be de- 
termined by the following equation 

N 

t t = t x + (6Kq - ^)(N - 1) - 4Kq £ n 

n -2 



< o\ co 



< .J- G\CO 


11 


-0 

where N is the nearest integral value of 2bX 7 + ^ an< ^ "*"°®e ^ I = 0 

0 

if ] | < 1. The quantity N represents the total number of trajectory 

segments which terminate on the switch curve prior to the motion which 
subsides along the switch curve. The last term in equation (12) repre- 
sents the time during which motion is along the switch curve and is the 
time to damp for a first-order linear system with a root equal to -l/K^. 

The fuel consumption is calculated from the summation of momentum 
or angular rate changes (eq. (8)) 


n' 

v t = e'o + 2N' (2aK* -6i) - ^aK 0 -£ n (13) 

n=i 

where N is the nearest integral value of 


-0i 

2aK 


- + 0.5 

6 


/ 

The quantity N represents the total number of trajectory segments 
with the angular note at the end of the segment of opposite sign from 
that at the beginning of the segment . For the remaining portion of the 
motion which is represented by a linear system, a sufficiently wide dead 
zone along the switch curve is assumed so that the fuel consumed is only 
that required to reduce the angular rate to zero. 


Minimum settling time system with nonlinear switch curve .- This 
system brings the vehicle to rest in a minimum time for an arbitrary set 
of initial conditions (ref. 8) as shown in sketch (d). 



Sketch (d) 
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The resulting motion obtained through use of this switch curve is given 
by the equation 


9 = -a sgnj 
The total settling time equation becomes 


{■ ♦ e -w) 


t. 


“t a 

The fuel consumption is simply 


+ 2 / ia + ^ 


2a' 


v t = at t = e 0 + 2J ae o +^f- 


( 14 ) 


( 15 ) 


Note that these equations represent the time and fuel required to bring 
the vehicle to zero error and error rate rather than to a given small 
error , as was the case for the previous system. However > the time and 
fuel difference is negligible for this system because of the large error 
range to be considered. 


Rate-limited system .- This system limits the maximum rate as shown 
in sketch (e) . 
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Sketch ( e ) 


The intersection at 0 2 ,0 2 was selected so that no overshoot occurred 
at the origin. A sufficiently wide dead zone was assumed so that coast- 
ing would occur during rate limiting. The equation for total settling 
time is 


, = 2a + !i + + 


2a Qr 


a 


(16) 
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The equation for fuel consumption is 

w t = e 0 + 2 e L (17) 

The rate limited system can be adjusted to yield a lower fuel con- 
sumption for a given settling time than other systems since fuel con- 
sumption is proportional to changes in q . This minimum is obtained 
by adjusting the rate limit to the minimum value possible to still 
achieve the desired settling time for the system. 

System with detector saturation .- The form of detector saturation 
being considered is shown in sketch (f). 



Sketch ( f ) 

Since the rate signal is assumed to be determined by differentiation of 
the detector signal, then the rate signal is unobtainable during detec- 
tor saturation. These combined effects may cause a considerable amount 
of overshoot, as shown in sketch (g), and result in a deterioration of 
settling time and fuel consumption . 
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Sketch ( g ) 


Point A is the last point the trajectory intersects on the 0 axis 
before the portion of the switch curve below the 6^ limit is inter- 
sected. Point B is the intersection of the saturated part of the 
switch curve with a selected switch curve for errors less than the de- 
tector saturation value. The time required to reach point A is 


t 


A 



M 

2 -7 2 /a D 

m=i 



+ It - (1 - 


(18) 


where 


M 


is the nearest integral value of 





Geometri- 


cally, M + 1 represents the number of intersections of the trajectory with 
the 0 axis . If a linear switch curve is used below saturation values , 
the total settling time is obtained by adding the times calculated from 
equations (l8) and (12) . The time from equation (l 1 )-) is used in place of 
equation (12) if a minimum time switch curve is used. The initial values 
to be used with either of the latter equations are the coordinates of 
point A. 
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The fuel required to reach point A is 

M 

= at^ = Q 0 + 2\^2a S 
m=a 

The total fuel expended to decrease to a given amplitude is obtained by 
the sum of equations ( 13 ) and ( 19 ) for the linear switch-curve case, or 
the sum of equations ( 15 ) and ( 19 ) for the minimum time switch curve. 

Timer system *- For cases in which detector saturation effects be- 
come important, a system which uses a timer to control torques can be 
desirable. Since changes in angular rates are proportional to time for 
a vehicle with only inertial moments, knowledge of the state of the 
vehicle can be obtained from a time and amplitude measurement. A timer 
can then be used to program switching times to bring the system to rest 

The motion for a system which would require only the direction of 
error beyond a given dead zone for the detector signal is shown by 
sketch (h). 




Sketch (h) 
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At point 1., a # constant positive torque is applied. After a period of 
time, t 2 = 2|0il/a, point 2 is reached. At point 2 the torque is switched 
to a negative value and is applied for (l/2 + */2/4)t 2 . At point 3; posi- 
tive torque is applied for Cs/2/4)t 2 , and the system is then at zero 
displacement and rate. For perfect time measurement and exactly equal 
positive and negative torques, the final switch point, 3j is the same as 
a point on the previous minimum settling time curve. The total settling 
time is 


Fuel consumption is given by the equation 

v t = at t = (3 + -Jz) |©1 

Since the system operates with an open loop for large portions of 
the trajectory, the effects of unequal magnitudes of the positive and 
negative applied accelerations, a x and a 2 , should be investigated. The 
resulting equations will also apply for cases with equal applied torques 
but with a constant disturbing torque. The equations used for calcu- 
lating settling time and fuel consumption will depend upon the relative 
magnitudes of the two accelerations . 

Motions of the system for one cycle with unsymmetrical acceler- 
ations are shown in sketch (i). 



( 20 ) 
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k 

9 

( 21 ) 6 


Sketch (i) 
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Note that the relative magnitudes of the applied accelerations cause the 
system to either overshoot (solid curve) or undershoot (dotted curve) 
the zero angle at the conclusion of the timed torque sequence. If 
undershoot occurs (point ^-a) , that is, if the system does not return to 
zero angle in the given time, negative acceleration is immediately 
applied until zero angle is reached. For either case, the cycle is 
then repeated until the residual angular rate at zero amplitude is within 
a desired value. In the subsequent equations to be given for this system, 
a small dead zone is assumed which represents the desired tracking 
accuracy, 0^. Since this value is very small in comparison with the 
initial motions, it is not shown in the sketches. 


Case I, overshoot (ax > a 2 ): For this case, the system will reach 
zero amplitude before the angular rate is zero . First consider only an 
initial angular rate 0 i • The settling time for the first cycle can be 
expressed in dimensionless form as 


at c 

"el 1 


a t 2 a tg a 1 4 

JeTT + TeTT + TeTT 


where 


( 22 ) 


at 2 




at 3 




at4 

FIT 




g ( e 3 ~ e d) a 

a 2 e x 2 



-l 


+ 



at 

3 




The ratio of initial to final rates is 



( 23 ) 
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The fuel consumption for the cycle is 

w c = 3l0il + 21 03 1 - l0 4 | 

The total settling time can he expressed as a geometric series of the 
settling times for each cycle so that the total torque-on time required 
to subside to a desired angular rate, 0 d > is 


at T _ at c /I - |ea/0i|\ 

IM ' Mil Vl - I MM (24) 

f 

1 

Note that the ratio B&Jq i does not represent, in general, a ratio that 
can be obtained from an integral number of cycles. However, if 0^ is 
small in comparison with the effect of the discrepancy in the final 

value of 6 on total time is negligible. The time required to cross 
the dead zone becomes significant as the angular rate becomes small. If 
^d and 0& are small in comparison with the maximum values of Q and 0 
obtained during the first cycle, then the effect of dead zone on settling 
time can be approximately expressed by the equation 


= th A - iMh h 

9 d A - I 04/01 1/ 


(25) 


For arbitrary initial conditions, the total settling time becomes 

tt = ti + tT + t z 
or 


where 


tt - ti + 



+ ^ A - I MM 

%/ v - I 04/01 1/ 



The fuel consumption is 


wt = 0 O + 


01 1 


+ w c 


1 - lMM\ 
.1 - I MAI/ 


( 26 ) 


(27) 
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Case II, undershoot (a x < a 2 ) : The first cycle is shown in the 

previous sketch. When undershoot results at the end of the last timed 
interval (point 4 a,) , opposite torque is applied immediately till zero 
amplitude is reached. The time for the first cycle is expressed as 


atr 


at? 


atr 


at 4a . at 4 

*r + 


10 : 


I 0 i I 


(28) 


where 


at 4a 
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1 0 i 1 

2a 2 


at 4 

__ 9 4a 
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1 ©I 1 

1 0i l 

a i 

040 , 

.a* 

n/X 

r-' 1 
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04 a B_'f 2(04a - 0d)a 2 


0 i a i 


a i 0 ; 


l - Si 

ap 


a ( 04 a - 0 d) a (03 - ®d) 


0 1 2 


0 i 2 


03 

01 


at 4 a 1 a 2 / a t 4 a \ 2 

| 0 x I 2 a \ 0 1 


The equations for the remaining quantities in equation ( 28 ) are the 
same as those for Case I. 


The ratio of final to initial rates is 


0 4 _ 9 40 , ^ 3 a at 4 

01 I 0 1 1 3 | 0 1 1 


(29) 


To find total settling time, equations (29) and (28) are substituted 
into equation (26) . The total fuel consumption is obtained by equation 

(27) with 

w c “ 3l$il + 2| 0 3 | + 2)040,1 + I 04 I (30) 
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Modified timer system .- As will be shewn in the discussion, the 
performance of the timer system deteriorates considerably if the applied 
torques are very unequal. One method for improving this limitation is 
to initiate the corrective torque sequence only when the error returns 
to the same sign that was present while the time measurement was made. 
This system is illustrated in sketch (j) . 



© 


Sketch (j) 

The operation of the system is as follows. The time to travel from 
point 1 to point 2 is measured with a constant positive torque 

applied. A negative torque is then applied to drive the system to point 
3. At point 3, a positive torque is applied for (l/2 +*j2/k)t s . At 
point 4, a positive torque is applied for ( v ?/4)t 8 . Again, differences in 
magnitude of positive or negative torques will cause the system to 
undershoot or overshoot the origin. 
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Case I, undershoot (a x > a 2 ) : For an initial angular rate the 

settling time for the first cycle can be expressed as 


cffcQ 0,^2 £L*t3 3.^4- ^^500 9,^5 

| 01 1 | 01 1 jfill I ©1 I I 01 I I 0l| 


( 31 ) 


where 
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The ratio of final to initial rates for the cycle is 


6 

0 i| 


* 50 , 
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Case II , overshoot (ai < a 2 ) : 
at c 


where 


e- 


at 2 at3 ai 4 at^ 

* I i" 1 : r I * 1 "i l * 1 
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a sj 2 
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us) 


2 2(e 4 + G d )a.‘ 


6 i 2 ai 


(32) 


( 33 ) 


The ratio of final to initial rates is 


0 1 


a l s ^5 
2 a | 6 2. | 


( 3 ^) 


The equations for the remaining quantities in equation ( 33 ) are the same 
as those for Case I. Total settling time end fuel consumption can be 
found in a manner similar to that used for the previous timer system. 

Rate system .- Another system that requires only a small error 
range for the detector is one that uses a lead network so that angular 
rate and angle can both be determined while the error is within the 
detector range. While a greater error range is required for this case 
than for the previous one in order to determine the rate, it still can 
be quite small. After the angular rate at zero error is determined, a 
timed torque sequence is initiated to brin£ the system to rest . A 
typical trajectory is shown in sketch (k) - 
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Sketch (k) 


The angular rate is measured at point 1. Positive torque is then 
applied for a time t 2 = (>.« - ^ , At point 2 , negative torque 


applied for a time t 2 = yl 
is applied for t 3 = 


The system is then at rest. The total 


settling time is 


t t = ( 1 + ^i) J|iL 


while fuel consumption is simply 


= at = (1 + n/ 2) 1 6i 
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Motions of the system for one cycle with unequal torques or rate 
measurement errors are shown in sketch (l). 



Sketch (l) 

A rate measurement error 6xxi> or torque errcr will cause the vehicle to 
again overshoot or undershoot zero angle at the conclusion of the timed 
torque sequence. The equations including these effects for settling 
time and fuel consumption will now be given. 

Case I, undershoot (a x > a 2 or \jhjn\ < 1 9 x | ) : The vehicle does not 

reach zero error in the computed time sequence so positive torque is immedi- 
ately applied to restore the error to zero. ?or an initial angular rate 
0 1} the settling time for the first cycle can be expressed in dimensionless 
form as 


at c at 2 a "^ 3 (x at 3 



( 37 ) 



where 


< o\co 
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where 


1 a x -J \ a i/ 


ms 

a i 


Fuel consumption is 


i- e,A 


w, = 0 n + 01 + w 


Cl- 




(4o) 


where 


w - ©i + @2 + 2 63 a +-^3 


Case II, overshoot (ag < a 2 or l_£jjg] > l^xl) : For an initial angular 
rate Qi, the settling time for the first cycle is 



(41) 


where 

• [f * \ 2 2'e s + e)a. 2 

at 3 _ a e 2 + — j i 

I ©i| a i|^i| V \ai ©i/ a l 0 l 2 

The expressions for the remaining quantities in equation (4l) are the 
same as for Case I. The ratio of initial to final angular rates per 
cycle is 


. • 

a i at^ 



Fuel consumption per cycle is 

v c =| e x | + 2\9 S \ - 1 e 3 | (42) 


Total fuel consumption and settling times are found by substituting 
equations (4l) and (42) into equations ( 39 ) ani (40), respectively. 


< -4 On CO 
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Limit-Cycle Characteristics 


■+ 

3 


Fuel consumption-- Once a system with adequate settling time has 
been determined, the limit-cycle characteristics must be investigated to 
establish that the fuel consumption rate will be sufficiently low. The 
fuel consumption per second can be expressed in terms of the average 
change in angular velocity by the following equation which is analogous 
to equation (8) 


W 

av 


Iw 

av 


I 1 


sp 


(*U) 


where w av is the average change in angular rate due to applied torque 
during limit-cycle motion. 


With no external torque present, the average change in angular rate 
is determined by the torque-on time, which will be determined in a later 
section, and the width of the coast zone, 0 CO ■ As will be shown sub- 
sequently, the limit cycle motion is symmetrical for a symmetrical 
linear switch curve system. In determining the fuel consumption rate, 
the total time will be approximated by the coasting time since the 
coasting time is much greater than the torque-on time. The average fuel 
consumption per second during limit-cycle motion can be expressed as 


= ( aT zh 

^co 


(MO 


Hence the torque-on time, , should be as small as possible, while the 
coast zone, 0 CO , should be as large as possible, consistent with track- 
ing accuracy requirements. 

In the presence of an external torque which is sufficiently large 
that a corrective pulse does not cause an overshoot of the dead zone on 
the opposite side, the fuel consumption is simply proportional to the 
acceleration produced by the external torque, 


This amount of fuel would have to be expended even for a momentum ex- 
change system, such as inertia wheels, in order to reduce the wheel 
speed periodically. For cases in which intermediate sized torques re- 
sult in both dead zone boundaries being intercepted, the fuel con- 
sumption depends more on the particular system, and relations for 
general initial conditions at the dead- zone boundary are more difficult 
to obtain. 
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Tor que- on time for system with line ar switch curve. - The limit - 
cycle characteristics of the linear switch curve system are of parti- 
cular interest since a linear switch curve could he used in conjunction 
with other systems in the low-amplitude limit- cycle range. The effects 
of hysteresis, time delays, time lags, and external torques on limit 
cycle characteristics will be investigated. The effects to be con- 
sidered are shown in sketch (m). 



Sketch (m) 

When a thrust- on signal is applied, a typical thrust transient response 
has the form indicated in sketch (n). 


Thrust 



— H f r |-» f » *"| f d \* 

Time 


Sketch (n) 

The thrust transient response has a rapid build-up and an exponential- 
type decay. The response can be approximated with sufficient accuracy 
by a step change in thrust and by the selection of time delays that keep 
the total area of the thrust-time curve the same. The time t r is the 
delay, after the signal-on instant, before the step thrust is applied. 
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The time t s represents the signal-on time after the step thrust is 
applied. The time tfl is the step thrust-on time representing the 
effect of thrust decay and thrust-off delay. 

Typical motion near the limit- cycle condition without external 
torque is shown in sketch (o). 



Sketch (o) 

Since only small values of torque-on time will he considered., the para- 
bolic torque- on trajectories approximate vertical straight lines. 
Although initially unsymmetric, the motion tends toward a symmetric 
limit cycle. It can be shown that the motion will converge to a 
symmetric limit cycle for sufficiently positive values of Kq. This 
convergence is due to small changes in torque-on time with initial 
values of 0 at the switching boundary. The equations for torque-on 
time attributable to system imperfections will be given subsequently. 

For the case with an external torque sufficiently large that only 
one switching boundary is reached, the motion becomes symmetrical only 
about the 0 axis, but the same limit-cycle torque-on time occurs. 

For cases with intermediate values of torque, the torque-on time varies 
but it will be close to that corresponding to the limit cycle motion. 
Thus an estimate of the torque-on time for the limit cycle without an 
external torque gives a good indication of limit-cycle characteristics 
with an arbitrary value of external torque. A further discussion of 
limit cycles in the presence of external torque will be given in a sub- 
sequent section. 

The coasting period is assumed sufficiently long for the signal 
from the differentiating network to reach a constant value equal to the 
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coasting angular velocity 0 O . For this initial condition, 0 O , the 
transient response of the output of the differentiating network, N(t), 
due to a negative step torque input is 


N(t) = Kq 


at 


+ aT 



e 



(46) 


In the equations that follow, it will he convenient to combine the time 
parameters with the angular acceleration, a. The signal-on time after 
the torque is initially applied, t s , (see sketch (n)) can he obtained 
from the following transcendental equation: 


ah + at-^o + at s 0 o 



( 47 ) 


The first two terms in equation (47) represent the change in control 
signal (sum of amplitude and differentiating network signals) required 
to initiate a cut-off signal. The remaining terns are the changes in 
the control signal while torque is applied. The total torque-on time, 
t, is obtained by adding the thrust decay time to the signal-on time. 


T “ t s + t d 


(48) 


Equations (4-7) and (48) are used to determine the torque-on time which 
occurs for an arbitrary initial angular rate, Q c , at the switch-curve 
boundary. The equations are valid for parameters for which torque-off 
occurs within the dead zone boundary. Two simplifications of equation 
( 47 ) are of interest. For the first case with relatively large hysteresis 
and a small time constant, the exponential term can be neglected. After 
making this simplification and combining equaticns (47) and (48) , the 
equation for torque-on time becomes 


P 

1 

c 

l 


aT - 0 O - oKq + at^ + 


(aK^ - 0 O ) 2 + 2ah + £at r 0 o + 


li 


(49) 


Results from this type of approximation were discussed in references 5 
and 6, For the second case, which is of interest for this report, a 
relatively large time constant is used since a large value of is 

desirable because of the low values of angular 1 ate encountered. The 
large time constant is also useful in reducing effects of noise. When 
the signal-on time is short in comparison with the time constant, the 
exponential term in equation ( 47 ) can be expand*:! into a power series* 

By including the first three terms of the series, the following expression 
for torque-on time is obtained from equations (^ 7 ) and (48) , 
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aT = 


0 O 2 + 2(2ah + at d 0 o ) (l + aJK^/aT) 
1 + aKo/aT 


+ at : 


(50) 


The torque-on time for the limit-cycle motion, ar^ , cam he obtained 
by use of the relation 

an = 20 o (51) 

After substituting equations ( 51 ) and (48) into equation (47) , a tran- 
scendental equation for torque-on time for the limit-cycle case is 
obtained. 


•i / . , .n , . 

— (^at r + at^; + ah - — - - aX^ 

Two approximate solutions of this equation can also be obtained. For 
the first case, with relatively large hysteresis and small time constants, 
equation ( 52 ) reduces to 

2ah - (at^) 2 + 2aK^(at^ + aT) ^ 

^ 2oKq - at r - at^ 


r l “ a ^d - a-T + aT e 


at,: 


aT z 

aT 


= 0 


( ‘SP ) 


For the present case with relatively large time constants, the expo- 
nential term in equation (52) is expanded into the first three terms of 
the power series and the resulting expression is 


aT z = at^ + (at r + at^) 


T 

2Kq 


+ 


at^ + (at r + at^) 



p 

(a.t d ) 




2Tah) 2 

*3 J 


(54) 


Equation (54) indicates that the thrust rise time has a smaller effect 
on T i than the thrust decay time since (t/K^)< < 1. If the thrust-rise 
time is approximated by the thrust decay time, equation (54) becomes 


aT Z = at d 



+ 


~(at d ) 2 T / 

f + 

L % ' 

V V 


2ahT 

- 


(55) 


When the external torque is sufficiently large that the switch 
curve on only one side of the dead zone is encountered, the motion 
approaches a limit cycle with symmetry only about the 0 axis. This 
motion will be discussed in more detail in a subsequent section. The 
previous equations for torque-on time can also be used in this case. 

The additional assumption is made that the effect of the small external 
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torque on the output of the differentiating network during the time, that 
the torque is applied is negligible. An additional equation of interest 
for the case with external torque is the one for determining the torque- 
on time during limit-cycle motion which will just intersect both switch- 
ing boundaries. The equation representing this limiting case of the 
intersection of the parabolic trajectory with the linear switch curve 
with a dead zone, 6 <± z > is 


a tj = - 2K^ a e + \j&e ® dz 


( 56 ) 


RESULTS AND DISCUSSION 


The predicted performance of the systems described in the previous 
section will be investigated. This will be done by a comparison of 
settling times and fuel consumption during recovery from the set of 
initial conditions described in sketch (a) in the previous section. In 
addition, limit-cycle characteristics (principally for the system with 
a linear switch curve) will also be discussed. Typical motions to- 
gether with resulting fuel consumption rates will be shown. In sub- 
sequent calculated results, it will be convenient to show variables in 
a normalized form, based on the desired control accuracy, 0 . The vari- 
ables given in the previously derived equations in dimensional form were 
normalized to a unit value of desired tracking angle by dividing them by 



Recovery From Initial Errors 


Linear switch-curve, minimum settling time, and rate-limited sys- 
tems .- The performance of the systems with linear and nonlinear (mini- 
mum time) switch curves, and the rate-limited system are shown in fig- 
ure 2. The linear and nonlinear switch curve systems will be discussed 
first. Typical phase plane motions for a specific level of acceleration 
are shown in figure 2(a). Similar curves were obtained for other ac- 
celeration levels and, for each acceleration, the value of K * for the 

linear switch curve was adjusted to give a minimum average settling time 
for the range of initial conditions shown in figure 2(a). In making 
this adjustment, the nonlinear switch curve for the same value of ac- 
celeration was used as a guide. The linear switch curve system requires 
a longer settling time than the nonlinear switch curve system, although 
the difference in fuel consumption is relatively small (figs. 2(c) and 
2(d)). The increase in time occurs principally during the latter part 
of the transient at low velocity levels and small average applied 
torques. At these low velocities, the linear system produces an expo- 
nential decrease in velocity with time with a decreasing average value 
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of torque, while the nonlinear system produces a linear decrease in 
velocity with full torque applied in one direction. For a settling 
time of about one minute, approximately 30 percent more fuel must be 
expended by the linear system than by the nonlinear system. The larger 
values of acceleration allow shorter settling times but require more 
fuel. Also, if very large accelerations were considered, fuel con- 
sumption during limit-cycle operation would become too great when 
effects of imperfections such as time delays were included. These 
limit-cycle characteristics will be discussed in a later section. 

As an example of the weight of the fuel required to recover from 
an initial error for a typical case, consider a vehicle with I = 

10 10 gm cm 2 , l = 75 cm, a desired tracking accuracy of 0 = 0.1 sec- 
ond of arc, and a relatively low efficiency fuel, such as water vapor, 

for which I = 50 seconds will be assumed. Note that for an acceler- 
sp 

ation of 36^/ sec 2 , a total thrust of only about 200 dynes is required. 
After a typical value of w t /0 d = 120/sec is selected, I is converted 

to cgs form, and these values are substituted into equation (8), a 
weight of 0.158 gm is obtained. If this correction is required once 
every other 100 minute orbit about each of two axes for one year, a 
total of 830 grams of fuel would be expended. 

Typical phase plane trajectories for the rate-limited system are 
shown in figure 2(b). For a given settling time, the fuel consumption 
can be considerably reduced by the use of rate limiting. The rate 
signal could be produced by either a lead network operating from the 
detector signal, or from a separate source such as a rate gyro if the 
rates are sufficiently large. 

An important effect to be considered is saturation of the detector 
since a large error range of 1000 to 1 is desired. The detector output 
is a signal assumed proportional to error angle. However, at some value 
of error input, the detector may saturate and its output will be limited 
at this value for increasing error angles. If a lead network is 
assumed for rate information, error rate is also only available when 
the detector is operating within its proportional limits. The effect 
of detector error saturation on the performance of the linear and non- 
linear switch curve systems is shown in figure 3* In the determination 
of the linear switch curves, the values of selected were those 

which gave minimum time response for the error and error rate range of 
interest below each error limit value. A large increase in settling 
time and fuel consumption occurs as detector saturation limits become 
lower. For the lower values of detector saturation, there is little 
difference in performance between the two systems since most of the time 
snd fuel is utilized in recovering to within the proportional range of 
the detector. With a detector limit of 6^/6^ = 50, a three minute 
settling time can still be achieved for the acceleration range shown. 
However, fuel consumption has increased by a factor of about five from 
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the unsaturated case to 4,300 grains. For a detector saturation of 
e L^d = the ^ ue l required for the system to settle within three 

minutes is reduced to about 2,300 grams. F.ence, a detector range of 
at least 50 6 ^ to 100 6 ^ would be required. For an application such 

as the orbiting astronomical observatory, this range may be difficult 
to obtain since problems of alinement dictate that the error signal be 
obtained from the primary optical system. 

Timer systems. - Two systems which can operate with a minimum 
detector range will now be discussed. In loth cases, a high resolution 
timer is used to apply a timed sequence of torques, open loop, between 
measuring intervals. The timer accuracies depend upon the residual l 

angular velocity desired, and could be as high as 10 milliseconds. The 
effects of torque magnitude errors will be discussed. As a limit cycle 
condition is approached, a second mode of operation generally must be 
used. Otherwise the requirements on system lags, etc., would become 
too stringent. 

The performance of the first timer system which requires only the 
direction of error from the detector is shewn in figures 4 and 5 - The 
system operates as shown in figure 4(a). Ihe time required to move 
from condition 1 ( 0 , 0 i) to condition 2 with positive torque applied is 
measured. Then opposite torque is applied for 85-3 percent of the first 
measured time. Then positive torque is applied again for 35*3 percent 
of the original measured time. If the applied torques in each direction 
are symmetrical, then final switching occurs as shown by the nonlinear 
switch curve in figure 2(a), and the vehicle is returned to rest. How- 
ever, if negative torque is greater than positive torque (ai>a 2 ), the 
trajectory tends to overshoot the origin ol the phase plane and the 
angular rate is too large when zero angle is reached. One cycle of the 
motion is completed. The cycle is then rejeated until the measured 
timed interval is small enough to indicate that the angular rate is 
reduced to a desired level. Note that quite accurate time resolution 
would be required as a limit cycle condition is approached. If negative 
torque is less than positive torque, undershoot occurs at the end of the 
35*3 percent time interval, and negative torque is then applied until 
zero angle is reached. The cycle is then lepeated. 

The ratio of final velocity to initial velocity for each cycle, 

84 / 81 , is a function of the torque error (fig. 4(b)). For torque errors 
greater than ± 7*9 percent from the average value for the ai> a 2 case, 
and torque errors greater than 11. 5 percenl for the ai< a 2 case, J$ 4 | 
becomes greater than |0*ij and the system is unstable. The variation 
of time per cycle with torque error is shovn in figure 4(c) and is seen 
to be fairly small. The effect of torque error on total settling time 
with = 0 is shown in figure 4(d). A large increase in 

settling time with increasing torque error is indicated. 


\ n 
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For a comparison with the previous systems, the performance of the 
timer system was computed for the previous set of initial conditions 
(fig. 5 ) through use of equations ( 26 ) and ( 27 ) with 0 ^/ 0 ^ = 0.25 and 

0^ = 1. The values shown for a given percent torque asymmetry represent 
the arithmetic average of the ai<a<a2 and ai>a>a 2 relative torque 
cases. When compared with figure 3 > the settling time and fuel con- 
sumption of the timer system without thrust error are seen to be 
approximately equivalent to the linear switch- curve system with detector 
limiting of about 150 0/0. With a ±2 percent thrust error, the per- 

formance is equivalent to the linear switch-curve system with about 
50 0/0 detector limiting. Hence the timer system becomes competitive 

CL 

only for low values of detector limiting. 

For large thrust errors it becomes more efficient to apply correc- 
tive thrust only when the error returns to the same sign that existed 
when the time measurement was made. This modified system is shown in 
figure 6 . A comparison with figure 4(d) indicates that while the time 
required by the modified timer system for the symmetrical thrust case 
is almost twice that of the timer system, the total time becomes equal 
or less when the thrust error becomes ± 1. 5 percent or greater. 

The results for the second timer system to be called the rate 
system are shown in figures 7 * 8 , and 9 with a typical trajectory 
shown in figure 7(a). The rate 6 1 is measured at 0 = 0 by use of a 
lead network. Positive torque is then applied for 1.706 9 i/a seconds. 
Then negative torque is applied for O .706 0 X / a second and the vehicle 
is returned to equilibrium. For symmetrical torques, final switching 
occurs at the conditions given by the nonlinear switch curve (fig. 2 ). 
Accurate time resolution is also needed for this case. For unsymmetri- 
cal torques, the vehicle is not returned to equilibrium but will tend 
to either overshoot or undershoot (fig. 7 (a)) in a manner similar to the 
timer system and the cycle must then be repeated. The rate system pro- 
duces a much smaller rate error per cycle ( 63 ) than the timer system 
for a given torque error. This difference results largely from a 
smaller error in the final switching instants (figs. 4(a) and 7 (a)). 

In addition to inaccuracies in torque magnitude, this system would also 
be subject to inaccuracies in rate measurement by the lead network 
(fig. 8 ). A comparison with the effects of torque error shows that a 
given percentage of rate measurement error causes a slightly greater 
increase in average settling time than the same percentage of torque 
error (figs. 7 (d) and 8 (d)). 

The performance of the rate system in recovering from the previous 
set of initial conditions (fig. 9 ) was calculated through use of equa- 
tions (39) and (40) with 0^/0^ = 0.25 and 0=1. The values shown for 

a given percent torque asymmetry represent the arithmetic average of the 
a i<a<a 2 and relative torque cases. A comparison with the linear 

switch curve system (fig. 3 ) indicates that without torque error, the 
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rate system performance was approximately equivalent to the linear 

switch-curve system with detector saturation of 6 /Q = 200. With a 

1/ d 

±2 percent torque error, the rate system performance was roughly equal 

to the linear switch curve system with detector limiting of 6 /6 = 100. 

1/ d 


Limit-Cycle Characteristics 


For any of the previously discussed systems, the limit-cycle 
characteristics must be investigated to see if a stable limit cycle 
exists and if fuel consumption is sufficiently low. The character- 
istics of a system with a lead network will be investigated primarily. 
The network could be used with any of the previously discussed systems 
since a proportional error signal would be available during motions 
close to limit-cycle conditions. 

Linear switch-curve system. - The influence of lead network param- 
eters, hysteresis effects, and time delays an limit cycle performance 
for an on-off system with linear switch curve will be investigated. 
Typical motions near the limit cycle with different amounts of external 
torque present are shown in figure 10. Witr no external torque 
(fig. 10(a)), the motion shifts toward a symmetrical limit cycle 
asymptotically. With an external torque sufficiently large so that only 
one switch curve is encountered (fig. 10(b)), the motion shifts toward 
a limit cycle which is symmetrical only about the 0 / 6 & axis. For 

intermediate values of torque (fig. 10(c))for which both switch curves 
are intersected, a true limit cycle is not attained, but the motion 
gradually hunts back and forth within a region determined by the switch 
curves. The time intervals during which torque is applied near a limit- 
cycle condition will be referred to as torque-on times. Since these 
times are very small, the parabolic torque-on trajectories approximate 
vertical straight lines. The variations in torque- on time near the 
limit cycle cycle are seen to be very small and they approach the same 
value for the limit cycles shown in figures (10a) and (lOb). 

The selection of torque-on time for the limit cycle motion depends 
upon the magnitude of external torque expected. As stated in the analy- 
sis section, fuel consumption for the motior shown in figure 10(b) 
depends only on the magnitude of the external torque. Hence it would be 
desirable to make the torque-on time as small as possible so that only 
one switch curve would be intercepted (fig. 10(b)) for a wide range of 
external torques, and therefore any fuel expended would be used only in 
counteracting the external torques. However, a smaller torque-on time 
for a limit cycle places more stringent requirements on system imper- 
fections, such as hysteresis and time lags. Typical variations in 
torque-on times for a limit cycle for cases with large time constants 
m the lead network (eq. (55)) are shown in figure 11. 
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As a compromise between these conflicting requirements, suppose for 
the example in figure 10 that a sufficiently small torque- on time is 
selected so that only one switch curve is encountered with external dis- 
turbances present as low as 0.005 6^/ sec 2 . The other system constants 

for the figure (10) example are = 5-0, 6 dz / e d = 1-5- Then, from 

equation ( 56 ), a normalized torque-on time, aTj , not greater than 
0.29 0 ^/sec should be used. From the previous discussion on recovery 

from initial errors, a typical value of applied acceleration is 

2-5 0^/sec 2 . Hence, the firing interval required is about 0.1 second. 

To achieve this pulse increment with the given value of acceleration, 
quite stringent requirements are placed on values of time lags and 
hysteresis. For instance, for K^/ T = 10, a combination of hysteresis 

of 0.08 9 ^ and a time delay of 0.032 second would be needed to achieve 

the desired pulse size. An additional effect to be considered in the 
selection of the dead zone is that, although only angular tracking 
accuracy is required, the linear switch curve will produce switching at 
different values of angular error, depending on the amount of angular 
rate present (fig. 10 ). This variation is particularly troublesome for 
motions with intermediate amounts of external torque present 
(fig. 10(c)). Hence, the switch curve dead zone should be selected so 
that 6 (sketch (o)) during limit cycle motion (without an external 

torque present) is somewhat less than the desired tracking accuracy, 6 ^ . 

This reduction in dead zone increases the fuel consumption for cases 
with zero or very small amounts of external torque. 

An estimate of the weight of the fuel consumed during the limit- 

cycle motion can be made for the previous example of I = 10 10 gm cm 2 , 

I = 50 seconds, Z = 75 cm, and 0 = 0.1 second of arc. If an accelera- 
S P d 

tion of 0.02 9 d per sec 2 due to an external torque is present, then the 

resulting fuel consumption (eqs. (45) and (* 1 - 3 )) is W &v = 2 . 63x10" 5 

gm/sec. The total fuel required for continuous two-axis control for a 
year for these assumed conditions would be 1650 grams. If no external 
torque is assumed, 0 CO = 0.09 second of arc, and aTj = 0.29 ©d then 

the resulting average fuel consumption (eqs. ( 44 ) and ( 43 )) is 

W av = 3-09X10 5 gm/sec. 


Incremental control system. - The performance of an alternate system 
for a case for which the angular rate becomes too small to obtain an 
adequate rate signal is shown in figure 12. In addition, this system 
can be used with dead zone values closer to the desired tracking angu- 
lar value since the tracking angle at which the pulse is fired does not 
depend upon angular rate. This incremental- type system uses pulses of 
fixed magnitude applied at discrete angle increments. To cope with the 
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effects of external torques, and still use a desired angle increment and 
pulse size, a memory of the two previous angle increments is required. 

As shown in the figure, a single pulse is applied at each angle 
increment that is of greater absolute magnitude than the previous one. 

If the same angle is reached again after a single pulse is applied, then 
a double pulse is used (fig* 12(b)). (For tne 0/0^ = 1 value, only 

single pulses are used). For both the time histories shown, the 
settling time was fairly long (about l-l/2 minutes). Note that if the 
incremental system is combined with the previously discussed on-off 
linear switch curve system, the limit-cycle torque- on time for the 
latter system must still be determined to insure that the resulting 
pulse is of the same order of magnitude as tne fixed pulse. 


CONCLUSIONS 


The performance of several systems using on-off reaction jets for 
precision attitude control of a satellite has been investigated. Both 
effects of recovery from large disturbances and of limit- cycle charac- 
teristics have been included. 

1. Detector saturation was found to have a major influence on 
settling time and fuel- consumption characteristics for the example con- 
sidered. For a lead network system, a detector range of at least 50 to 
100 times desired tracking accuracy was required to give satisfactory 
values of settling time and fuel consumption. 

2. Rate limiting was an attractive way to reduce fuel consumption 
and still obtain a specified settling time. When a lead network was 
used, no signal was available during detectcr saturation. Hence an 
auxiliary source of rate signal for limit detection would be needed. 

3. The performance of a system designel for use with a minimum 
detector range by means of a timed sequence of torques was also investi- 
gated. This system required a high resolution timer of up to 10 
milliseconds accuracy. Although this system required a longer settling 
time and higher fuel consumption rate than tne lead network system with- 
out detector saturation present, it achieved better performance than the 
lead network system with low detector saturation values. However, dif- 
ferences in positive and negative applied torque levels were found to 
have a large effect on settling time and fuel consumption of the timer 
system; for example, for a torque difference of ±2 percent, settling 
time was doubled. 

h. Another system which used a combination of a timer and a rate 
measurement only during motion over a small detector range achieved 
settling times about one half of those for tne timer system. 
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5 . Fuel consumption was low during limit-cycle operation of a 
system with a linear switch curve, and realistic values of hysteresis, 
lead network constants, dead zone, and thrust time delays, provided the 
applied torque level was moderate. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif,, May 9, 1961 
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APPENDIX 


DERIVATION OF EQUATIONS FOR SETTLING 
TIME AND FUEL CONSUMPTION FOR 
LINEAR SWITCH- CURVE SYSIEM 


As stated in the analysis section, the motion for this system 
(sketch (c)) may be separated into three parts: 

(1) The time required to reach the switch curve initially from a 
given set of initial conditions, 

(2) The time required during which the trajectory overshoots the 
switch curve when it is intersected, and 

(3) The time during which the motion subsides along the switch 
curve until the desired attitude angle is reached. 


The time for the first part of the motion can be determined through 
use of the equations for the switch curve and the trajectory segment 


ti 




(Al) 


The angular rate at the first intersection of the switch curve is 


at-i 


(A2) 


For the second part of the motion, there are a series of trajectory 
segments which start and terminate at the switch curve. The time 
required for the nth trajectory segment can be expressed in terms of 
the angular rate at the beginning of the trajectory segment as 


t 


n 


+ 


20 


n-i 


(A3) 


with n > 2. The equation for 0^ at the enc. of the nth trajectory 
segment is 


e n I = | - 01 - 2aK 0 *(n-l) 


(A4) 


with n >1. For a given value of 6 1 which corresponds to the first 
intersection of the trajectory with the switch curve, the total time for 
motions represented by n-1 trajectory segments can be obtained by 
combining equations (A3) and (a4) 


< -=t o\co 
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The total number of segments, N, is determined by the requirement that 
the motion be reduced to small enough values so accelerations commanded 
by the switch curve will be within the acceleration capability of the 
vehicle. Geometrically, this condition occurs when the magnitude of the 
slope of the phase-plane trajectory is equal to or greater than the 
magnitude of the slope of the switch curve. The maximum angular rate 
for which this occurs is 



(A6) 


If equations (A6) and (AU) are combined, N is seen to be the integral 
value of n for which the following relation holds 

I -6>! - 2aKg(n-l) | < aK^ 

or N is the nearest integral value of (-0i/2aK^) + 1. 


Part three of the motion is the average motion that is commanded 
by the switch curve. A sufficiently wide dead zone along the switch 
curve is assumed so that corrective torque is applied intermittently in 
only one direction to reduce the angular rate to zero. Hence, the time- 
averaged value of torque is less than the steady torque- on value. The 
resulting average motion is described by a linear first-order differ- 
ential equation with a root equal to -1/Kq. After equation (At) is 

used to determine the initial angular rate for this motion, 0^, and with 

the given desired tracking angle, 0^, the third portion of the settling 

time, t , is determined by the equation 
b 


t 


L 


K e log e 



K e e i 


2aK- 2 (N - 

U 



(A7) 


where log e | | = 0 if | | <1. The total settling time is given by 

the sum of equations (Al),(A5), and (A7). 


= + ( 6K e - fh - i: 


N 


-4K • E n + K ’log 


n=2 


1_ 

0 _ 


[- 


K-0! - 2aK- 2 (W - 1 


)]| (A8) 
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In determining fuel consumption, the total number of switchings 
must be defined differently than for the settling time case since fuel 
consumption depends only on the summation of maximum velocities 
obtained. As shown by the trajectory which starts at 0 O *, ©o* in 
sketch (c), the last two switchings may occur with the same sign of © 
and hence the final switching should be excluded for this type of tra- 
jectory. This effect can be included through use of equation (A4) and 
by defining the final trajectory segment N' as the minimum value of 
n for which 


-01 - 2aK-(n - l)<2aK- 

V V 

A 

or N is the nearest integral value of ( -©/2aK •) + 0-5. The normalized 4 

U 

fuel weight is calculated from the sum of the aigular rate changes. 


IT 

w t = ®o + 2 £ |®n 
n=i 


(A9) 


After the substitution of equations (a 4) into equation (A9), the 
following expression is obtained. 

W ' 

w^_ = 0 O + 2N ( 2aK ^ - 0 X ) - 4 .iK^ ^ n (A10) 

n=i 



< -d- O CO 
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Iritial conditions 


^Nonlinear 
0 + 8.0 0=0 


(a) Typical phase plane trajectories for minimum time linear and nonlinear 

switch curve systems* 


e/e 6 0 


-80 ' — 
-200 




200 400 600 800 1000 


(b) Typical phase plane trajectories for rate-limited system. 

Figure 2.- Comparison of minimum settling time linear and nonlinear 
switch curve systems , and rate-limited system. 












(b) Typical phase plane trajectory for amplitude limit of 0 L /0 d = 50. 

Figure 3.- Effect of detector saturation on minimum settling time linear 

and nonlinear switch curve systems. 
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(a) Normalized phase plane trajectory; torque error = ±2 percent. 



Figure 4,- Effect of torqpie error on timer system response. 














(a) Typical phase plane trajectory.; ai/a 2 = I.02/0.98. 



Figure 5 »- Effect of torque error on timer system response to initial 

conditions. 
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(a) Normalized phase plane trajectory; torque error = ±2 percent. 



Torque error, percent a 

(h) Total settling time for an initial rate error. 
Figure 6.- Effect of torque error on modified timer system. 


00 VO 4 =- > 
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(c) Time, per cycle. 



0 ±1 ±2 ±3 ±4 ±5 ±6 ±7 

Torque error, percent a 

(d) Total settling time for an initial rate error. 


Figure 7.- Concluded. 
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0 ±1 ±2 ±3 ±4 ±5 ±6 ± 


Measured rate error, percent 


(b) Velocity decrement, per cycle. 

Figure 8 .- Effect of rate measurement error on rate system response. 







la r acceleration, a/0<j, per sec 

elative fuel consumption. 


9. - Concluded 
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(c) External disturbance; a e / 0 ^ s= 


Figure 10 . - Cone luded . 





















